'''
Created on Jun 1, 2009

@author: morales
'''

import msTools
import epochs

NUM_OF_SIMULATIONS = 1
NUM_OF_CHROMOSOMES = 10

N = 1e4
MU = 1e-8
LENGTH = 25000
THETA = 4.0 * N * MU * LENGTH

N1 = 1e3
T1 = 5000
myCommand = '%d %d -t %f ' % (NUM_OF_CHROMOSOMES, NUM_OF_SIMULATIONS, THETA)

myEpoch = epochs.Epoch2_CC( T1, N1, N )
myCommand += myEpoch.getParameterString()

print msTools.runMS( myCommand )
